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ABSTRACT 

Context. Currently, the thermal emission from exoplanets can be measured with either direct imaging or secondary eclipse mea- 
surements of transiting exoplanets. Most of these measurements are taken at near-infrared wavelengths, where the thermal emission 
of these planets peaks. Cool brown dwarfs, covering a similar temperature range, are also mostly characterised using near-infrared 
spectra. 

Aims. We aim to show how thermal radiation in brown dwarf and exoplanet atmospheres can be scattered by clouds and haze and to 
investigate how the thermal emission spectrum is changed when different assumptions in the radiative transfer modelling are made. 
Methods. We calculate near-infrared thermal emission spectra using a doubling-adding radiative transfer code, which includes scat- 
tering by clouds and haze. Initial temperature profiles and cloud optical depths are taken from the drift-phoenix brown dwarf model. 
Results. As is well known, cloud particles change the spectrum compared to the same atmosphere with the clouds ignored. The clouds 
reduce fluxes in the near-infrared spectrum and make it redder than for the clear sky case. We also confirm that not including scatter- 
ing in the spectral calculations can result in errors on the spectra of many tens of percent, both in magnitude and in variations with 
wavelength. This is especially apparent for particles that are larger than the wavelength and only have little iron in them. Scattering 
particles will show deeper absorption features than absorbing (e.g. iron) particles and scattering and particle size will also affect the 
calculated infrared colours. Large particles also tend to be strongly forward-scattering, and we show that assuming isotropic scattering 
in this case also leads to very large errors in the spectrum. Thus, care must be taken in the choice of radiative transfer method for 
heat balance or spectral calculations when clouds are present in the atmosphere. Besides the choice of radiative transfer method, the 
type of particles that are predicted by models will change conclusions about e.g. infrared colours and trace gas abundances. As a 
result, knowledge of the scattering properties of the clouds is essential when deriving temperature profiles or gas abundances from 
direct infrared observations of exoplanets or brown dwarfs and from secondary eclipse measurements of transiting exoplanets, since 
scattering clouds will change the depth of gas absorption features, among other things. Thus, ignoring the presence of clouds can 
yield retrieved properties that differ significantly from the real atmospheric properties. 

Key words. Radiative transfer - Methods: numerical - Planets and satellites: atmospheres - Stars: low-mass, brown dwarfs - Infrared: 
planetary systems 



1. Introduction emission spectrum of the planet compared to the case without 

scattering. This scattering behavi our is evid ent for therma l emis- 

In recent years characterisating the atmospheres of exo- sion from, e. g.. Venus (Grinspoon et al. 199 3: Tsang et al.^OOSl) 

planets has started to become possible using spectroscopy, and Jupiter (ICarlson et al. 1991 and can be expected to be pro- 

The planets for which measurements at multiple wavelengths nounced at the short infrared wavelengths where hot exoplanets 

are available now are either large transiting planets (e.g. emit most of their thermal radiation. 
iDeming & Seager 2009 , and reference s therein ), or large planets 



far from their star (e.g. iPatience et al . 2010: B owler et al. 2010t Clouds have alr eady been detecte d on brown dwarfs 
iJanson et al. 2010| ). For these planets the thermal emission can (Ts uii et al. 1996t iLeggett et al. 19981: lAUard et al. 2001 ' 



be measured, either from its secondary eclipse in the former Ackerman & Marlev 200 It ICushing et al. 200a 



case, or by direct imaging in the latter case. These planets are ei- Burgasser 2003' and have been suggeste d for young plan 

ther very close to their star (transit measurements) or very young ets orbiting at large orbital distances (e g. iBowler et al. 20101: 

(direct detection), meaning that they tend to be hot and their Janson et al. 2010: Bonnefov et al. 20101) . Also, transit tran- 

thermal emission peaks at near-infrared wavelengths. Therefore, mission measure ments of hot exoplan ets have revealed 

most measurements of the thermal emission, from direct imag- hazes (iLecavelier Pes Etangs et al. 20081 : iPont et al. 20081: 



ing or secondary eclipse, are taken at near-infrared wavelengths. ISing et al. 20091 IWood et al. 2011h . " However, compar 



Along with reflecting starlight, infrared light from the planet it- isons of atmospheric models with infrared secondary 
self can be scattered when clouds or hazes are present in the eclipse measurements that probe the dayside the rmal emis 



planet's atmosphere. This scattering of the light can significantly sion of transiting exoplanet s (e.g. [B urrows et al. 2008a ; 



alter the heat balance in the atmosphere and change the thermal iKnutson et al. 200"9t iMadhustidhan & Seager 2009 ; 
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ISwain et al. 20091: ICroU et al. 201 Oh have so far neglected 
the possible presence of clouds or haze. 

At the relatively cool temperatures of substellar atmo- 
spheres, dust and cloud formation is predicted to significantly 
affect the thermal structure of the atmosphere. Some numerical 
algorithms are available that predict cloud properties in brown 
dwarfs, given a certain set of input parameters, such as the ef- 
fective temperature and the gravity. These models probe temper- 
atures similar to hot giant exoplanets. The available models can 
be divided i nto the ones that use a ki netic appr oach for dust for- 
mation (e.g. lWoitke & Helhng 2003r2 004; He lling et al. 2008ct) 
and the ones that assume phase equilibrium (e.g . 
Tsuii et al. 1996t lAckerman & Marlev 20011: lAllard et al. 20011: 



Barman et al. 200111 20051 iBurrows et al. 20061) . A recent 
comparison shows that the resulting atmospheres can differ 
substantially between the different brown dwarf cloud models , 
given the same set of input parameters (Hell ing et al. 2008al) . 
although the models that include both chemistry and opacity in 
the cloud formation show a relatively good agreement (F. AUard, 
pers. comm.). 

The radiative transfer methods that are used to calculate 
spectra (SEDs) and heat balance in the available brown dwarf 
models also differ. The equation of radiati ve transfer in a plane- 
parallel atmos phere can be written as (e.g. lChandrasekhar 19601; 
iHubenv 2003h : 



div . „ 

H— ^ ly-Sy 

dTy 



(1) 



Here, ji is the directional cosine, ly is the wavelength- 
dependent radiance, Ty is the optical depth, and 5v is the source 
function. In a planetary atmosphere all these normally vary with 
altitude or pressure and Eq. [T]has to be solved taking into account 
the entire atmosphere to yield the radiation that is emitted into 
space. The source function generally has contributions from both 
thermal emission of the local atmosphere itself and from scat- 
tering of radiation coming from other parts of the atmosphere. 
In local thermal equilibrium, the emission source function is 
the Planck function for black-body radiation. By. The scattering 
source function can be written as: 



PylydQ. 



(2) 



In the scattering source function the radiation ly coming in 
from different directions are integrated over solid angle Q and 
are weighted according to the phase function Py. The two contri- 
butions to the total source function are then weighted according 
to the single-scattering albedo 



5v = (1 - m)By + W()5v,I 



(3) 



The single-scattering albedo is the fraction of the light that 
is scattered when initially reaching a particle or layer in the at- 
mosphere. The remaining fraction is absorbed. The term single- 
scattering albedo should not be confused with 'single-scattering' 
methods of radiative transfer, where only one order of scattering 
is taken into account when solving Eq. 1. In the method used 
here, and also in the literature referenced here, multiple scat- 
tering radiative transfer is used, i.e. light is scattered more than 
once between different particles or atmospheric layers. In stellar 
literature the collisional destruction probability, e = 1 - wq, is 
commonly used instead of the single-scattering albedo. 

The equation of radiative transfer is solved differ- 
ently in different brown dwarf models. The phoenix 



model dHauschildt 1992t I AUard & Hauschildt 19951: 
iHauschildt & Baron 20061) increases efficiency of the cal- 
culations by assuming isotropic scattering, i.e. the phase 
function in Eq. |2] distributes the light equally in all directions 
(P. Hauschildt, pers. comm.). The equation of radiative transfer 
is then solved iteratively at multiple angles through the atmo- 
sph ere, using a method base d on approximate A-iteration (e.g. 
Hu benv 20031) . iBurrows et a l. (2006^ 20 08a) use a version of 
the TLUSTY code ( Hubenv & Lanz 1995h to iteratively compute 
the radiative transfer within the atmosphere. There radiation is 
calculated at multiple angles, with the phase function param- 
eterised in terms of the asymmetry parameter g (A. Burrows, 
pers. comm.). A value of g - \ would indicate perfectly 
forward-s cattering particles, whereas^ - indicates isotropic 
scattering) Ac kerman & Marlev (20011) use a two-step approach 
in calculating the spectrum and radiation field. Initially, they 
calculate the internal radiation field for the source function {ly 
in Eq. |2| at eight different angles in each atmospheric layer 
For the scatterin g term only they use the two-stream method of 
iToon et al. (19891) (M. Marley, pers. comm.) as the background 
source from which the scattering from particulates is computed. 
The two-stream approximation gives an exact solution for the 
case where radiation only consists of an up-going stream and a 
down-going stream. In this case the phase function can again 
be described using g, at the cost of a drop in accuracy. The 
two-stream source function method itself generally g ives error 
less than -10-20% for a single layer (iToonetal. 1989R . With the 
source function fully described, Eq. 1 can simply be integrated 
numerically to yield to radiation at the top of the atmosphere 
at multiple angles. Because Ackerman and Marley only use 
the two-stream method to compute the scattering source term 
and not for calculating the final spectrum and radiation field, 
their method is presu mably somewhat more accurate than the 
two-stream method of lToon et al. (19891) . 

Although scattering of the thermal emission is included in 
the previously discussed models, the various assumptions or ap- 
proximations that are currently used to model the scattering may 
lead to significant errors in the calculated spectra and hence in er- 
rors on atmospheric properties derived from observations. Here, 
we will calculate near-infrared thermal emission spectra for spe- 
cific model brown dwarf atmospheres, which are similar in tem- 
perature regime as hot exoplanets. Our calculations accurately 
include single and multiple scattering of the thermal radiation 
using scattering particles with angularly resolved phase func- 
tions. Hence, we will make an assessment of the effect of scat- 
tering by clouds and haze on the thermal emission spectra of hot 
exoplanet atmospheres and can compare spectra with different 
scattering assumptions in isolation, instead of comparing spec- 
tra from different cloud models that couple physics, chemistry 
and radiative transfer. 



2. Model atmospheres 

In this paper, we do not self-consistently calculate tempera- 
ture and composition profiles of the model atmospheres with 
our radiative transfer code, but instead use model atmospheres 
calculated with the self-consistent drift-phoenix code. With 
the DRIFT-PHOENIX code, atmospheric cloud particles are formed 
when rising trace gases cool down and seeds are formed. 
The particles are transported through the atmosphere, accu- 
mulating more material and/or evaporating again. The result- 
ing cloud particles usually consist of mixtures of several solids 
and are therefore called 'dirty' particles. Details about the ac- 
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tual computations perf o rmed by drift-pho enix can be found in 
iHeUing et al. (2008blch : rWitte et al. (20091) . 

In this paper, we use two model atmospheres that are cal- 
culated using DRIFT-PHOENIX, assuming solar metallicity and a 
gravitational acceleration of 10 ms"^. The model atmospheres 
only differ in the given effective temperature (Teff=1500 K and 
2000 K, respectively). The calculated temperatures profiles and 
cloud particle sizes and number densities for the two model at- 
mospheres are shown in Fig. [T] We use the vertical profiles of 
the temperature, trace gas mixing ratios (defined as the trace gas 
number densities divided by the total gas number density), cloud 
particle number densities and microphysical properties (particle 
size and composition) versus pressure as input for our radiative 
transfer calculations. The size distributio ns of the particles i n 
DRiFT-PHOENix are described by moments (iHelling et al. 2008cl) . 
but here we will only use the mean particle size as a function of 
pressure from drift-phoenix. 




0.001 0.010 0.100 1. 0001 0"'°1 0° 1 ° 1 " 1 0"M 0° lO'^ 500 1000 1500 2000 2500 3000 
Particle size (micron) Number density (cm"^) Temperaure (K) 

Fig. 1. Mean particle sizes, number densities and temperature 
profiles for our two model atmospheres. In the atmosphere with 
Jeff - 1500 K (reif = 2000 K), particle condensation starts at a 
pressure of 10"^ bar (2- 10"^ bar). The horizontal lines in the left 
panel indicate the boundaries of the four particle types that we 
define for our radiative transfer calculations (see text and Fig. |2]), 
which are different for the two atmospheres. 

The vertical profiles provided by drift-phoenix consist of 
~400 homogeneous layers. To reduce radiative transfer compu- 
tation times, we reduce the number of homogeneous layers to 
40 by interpolating the temperatures, pressures and gas abun- 
dances, and binning the optical thicknesses of the clouds. This 
significantly reduces the number of layers, while still resolv- 
ing the different cloud layers in the atmosphere. In addition, 
DRiFT-PHOENix predicts particles that have compositions and sizes 
that that are altitude dependent. As a result, the absorption and 
scattering properties (single scattering albedo, extinction cross- 
section and scattering matrix) of the particles depend not only 
on the wavelength but also on the altitude. Taking these vari- 
ations fully into account would require wavelength dependent 
calculations of the optical properties of 80 different particles (40 
layers and two model atmospheres). Instead, we define four dif- 
ferent particles for each model atmosphere (thus, eight in total), 
corresponding to four altitude ranges in each model atmosphere. 
Within each altitude range, the particle type is the same, only the 
particles' number density varies between the layers. Considering 





1500 K 


2000 K 


High silicate haze 


0.05 


0.03 


Low silicate haze 


0.07 


0.10 


Iron cloud 


0.22 


0.07 


AI2O3 cloud 


0.26 


0.25 



Table 1. Assumed particle sizes (jum) for the four different par- 
ticle types in the two model atmospheres. 



High Si haze Low Si Haze Iron cloud AI2O3 cloud 



Ti02 


0(0) 


0(0) 


1(2) 


4(4) 


Mg2Si04 


24 (23) 


30 (31) 


%(16) 


0(1) 


Si02 


22 (23) 


19(18) 


3(4) 


0(1) 


Fe 


15 (15) 


16(18) 


53 (35) 


4(10) 


AI2O3 


3(3) 


3(3) 


20 (29) 


92 (82) 


MgO 


1(8) 


6(6) 


W(8) 


0(2) 


MgSiOs 


29 (28) 


26 (24) 


5(6) 


0(0) 



Table 2. Assumed volume fraction (%) of constituents of the 
dirty particles for Teft- = 1500 (2000) K. 



the almost total lack of constraints on clouds on exoplanets from 
current observations and condensation models, this reduction of 
details seems reasonable for this exploratory study that focuses 
on the radiative transfer in a given model atmosphere. 

Our different particle types are based on the changes of par- 
ticle composition in the atmosphere. All particles produces by 
DRIFT-PHOENIX are altitude dependent mixtures of various solids. 
Figure |2] shows how the composition of the particles varies 
across the two model atmospheres. As can be seen, across large 
regions of each atmosphere, the particle composition only varies 
slightly with pressure. At a few altitudes, however, the particle 
composition changes rapidly. We choose these altitudes as the 
boundaries between our particle types. Based on these bound- 
aries, the particle types in the two model atmospheres have the 
same composition but different sizes (see below). We distinguish 
the following four particle types (from the top of the atmosphere 
to the bottom), named after their main constituent; a high silicate 
haze, a low silicate haze, an iron cloud, and an aluminium oxide 
(AI2O3) cloud. 

Although particle composition stays roughly constant with 
altitude across the four regions in each model atmosphere, par- 
ticle size does not (see Fig. [1). Since particle size is very im- 
portant in determining the scattering properties of particles, as- 
suming an altitude independent particle size across a region 
might have some consequences for our radiative transfer re- 
sults. Fortunately, the particles that drift-phoenix produces tend 
to be very small, especially in the upper atmosphere (above ~ 1 
mbar), where particle size changes most rapidly. When particles 
are much smaller than the wavelength, however, their extinction 
cross-section scales with the cube of their radius, wh ile its spec- 
tral dependence does not change (iHanel et al. 2003h . Hence, for 
the lower and upper silicate haze, we take into account the alti- 
tude dependent particle size in the calculation of a layer's op- 
tical thickness, but leave the spectral dependence unchanged. 
Lower in the atmosphere, the particles are larger, but their size 
also changes less with altitude. Hence, for the atmospheric lay- 
ers containing the iron and AI2O3 clouds, we assume a constant 
(the mean) particle radius. The assumed particle sizes used in the 
Mie calculations are given in Table 1. These coincide with the 
average of the mean particle size over the pressure range consid- 
ered. Table 2 summarises the composition of the dirty particles 
for each of the particle types. 
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1500 K 2000 K 




Volume fraction Volume fraction 

Fig. 2. The volume fraction of the different constituents of the 
dirty particles produced by drift-phoenix for the two model at- 
mospheres. The horizontal dotted lines indicate altitudes where 
the composition of the particles changes rapidly, which is where 
we have chosen the boundaries between the four different parti- 
cle types. The particle types can be described as: 1 - high sili- 
cate haze; 2 - low silicate haze; 3 - iron cloud; 4 - AI2O3 cloud. 
The altitude independent composition that we assume for each 
particle type is indicated by vertical dashed, coloured lines (see 
Table 2). 



3. Radiative transfer algorithm 

Using the model atmospheres produced by drift-phoenix 
as input, we calculate disc-integrated (spatially unresolved) 
near-infrared flux spectra with a doubling-adding algorithm. 
The doubling-adding method is often used in calculations 
of sunlight that is reflected by and/or transmitted th rough 
planetary atmospher es, see e.g. [Hansen & Travis (1974 ) and 
Ide Haan et al. (1987F) . It allows the accurate solving of the equa- 
tion of radiative transfer along paths in a model atmosphere that 
consists of a stack of plane-parallel, homogeneous layers con- 
taining scattering and/or absorbing gaseous molecules and/or 
aerosol particles. Vertically inhomogeneous model atmospheres 
are created by stacking different homogeneous layers. As de- 
scribed in Sect. 2. 1, our model atmospheres consist of 40 layers. 

As described in detail by e.g. Ide Haan et al. (19871) . the 
doubling-adding method makes use of the fact that if one knows 
the reflection and transmission properties of two adjacent atmo- 
spheric layers, one can straightforwardly calculate these proper- 
ties for the combined layer For each atmospheric layer, the cal- 
culation of its reflection and transmission properties starts with 
calculating them analytically for a horizontal slice of the layer 
with a very small optical t hickness. Then, u sing the so-called 
'doubling' -equations (see d e Haan et al. 1987h . the properties of 
a layer with twice the inital optical thickness are calculated. This 
is repeated until the required optical thickness is reached. The re- 
flection and transmission properties of the whole atmosphere are 
calculated by combining the pro perties of the indivi dual layers 
using the 'adding' -equations (see lde Haan et al. 1987i) . 

We use a special version of the doubling-adding algo- 
rithm, namely the so-called 'internal sources' algorithm of 
IWauben et al. (1994|), which fully includes multiple scattering 
of thermal radiation emitted by the planet's atmosphere. In this 
version, the temperature of each atmospheric layer (see Fig. 1) 



determines the amount of thermal radiation that is emitted in the 
layer The emission itself is isotropic. However, scattering of this 
radiation within the atmospheric layer and/or within other lay- 
ers will influence the angular distribution of the radiation as it 
emerges at the top of the atmosphere. We calculate the emerging 
thermal radiation along twenty angles between 0° (towards the 
zenith) and 90° (parallel to the atmosphere). Note that since our 
atmospheric layers are horizontally homogeneous, the emerging 
thermal radiation is independent of the azimuth angle. Next, as- 
suming the model planet is locally plane-parallel (thus that the 
mean-free path of the photons is small enough to ignore atmo- 
spheric curvature), we calculate disc-averaged thermal emission 
spectra by integrating the locally emitted spectra over the disc as 
follows 

</v) = 2 f ^ihd^ (4) 

Jo 

The cosine of the emission angle (yu) is included in the inte- 
gration to account for the spherical shape of the planet. 

To calculate the scattering and absorption of radiation with 
wavelength A within each atmospheric layer, the doubling- 
adding algorithm needs to know for each layer and at the given 
A: its extinction optical thickness (the sum of the scattering and 
absorption optical thicknesses), and the single-scattering albedo 
(the ratio of the scattering optical thickness to the extinction op- 
tical thickness) and the scattering phase functiorQof the mixture 
of particles and gas molecules in the layer 

Our model atmospheres not only contain the drift-phoenix 
particles, but also gas molecules. The altitude variation of the 
different atmospheric gases is produced by the drift-phoenix 
code. The absorption coefficients for H2 O and CO wer e calcu- 
lated using the HITEMP database ( Rothm an et al. 20Toh assum- 
ing Voigt line shapes. Other molecules that have strong spectral 
features in the near-infrared, such as CO2 and CH4 are predicted 
by DRIFT-PHOENIX to be present in only very small volume mix- 
ing ratios (generally less than 10"^ for CO2 and several orders of 
magnitude less for CH4), and should affect the calculated emis- 
sion spectra only marginally. Because the focus of this paper is 
on the scattering effects of clouds, and not so much on identify- 
ing gas signatures or comparisons with observations, we ignored 
the absorption coefficients of these low-concentration gases. 

To efficiently include the gaseous absorption in our calcu- 
lations of emission s pectra, we used the correlated-A: method 
(iLacis & Oinas 199T1) . We first calculated the absorption coef- 
ficients of the main gases at a high spectral resolution at a grid 
of 20 temperatures between 500-4000 K and 20 pressures be- 
tween 10"^ - 100 bar, which cover the relevant temperature- 
pressure space for hot exoplanetary atmospheres. From these 
high-spectral resolution spectra, we calculated and tabulated 
correlated A:-distribution coefficients that we interpolated to the 
temperature and pressure of each atmospheric layer Also in- 
cluded in our calculations is collision-ind uced absorption by Hg - 
H2, using the a bsorption coefficients of iBorvsow et al. (200 Ih ; 
iBorvsow (2002h . as well as Rayleigh scatte ring by H 2-H2, for 
which we use a depolarisation factor of 0.02 (llrwin 2 009). 

We calculated the wavelength dependent refractive indices 
of the various types of dirty particles that occur in our model at- 
mospheres from the refra ctive indices of the ir constituents using 
effective medium theory (iBosch et al. 2000i) . Then, for each type 
of particle, we calculated its extinction cross-section (i.e. the sur- 
face area with which one particle affects the light), its single- 

' In case polarisation is in cluded, the 4x4 phase matrix is required, 
see e.g. Ide Haan et al. ("1987h 
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scattering albedo and ph ase function across the required spectral 
range using Mie theory (Ide Rooij & van der Stap 19841) . assum- 
ing spherical particles with radii as given in Table 1 . Calculated 
extinction cross-sections and single-scattering albedos of the 
particle types occuring in our two model atmospheres are shown 
in Fig.[3]for the 1500 K atmosphere and in Fig|4]for the 2000 K 
atmosphere. As can be seen in these figures, none of the con- 
stituents of the dirty particles show any noticable absorption 
features below 5 jum, which results in extinction cross-sections 
that vary smoothly with the wavelength. The AI2O3 particles, 
although only micron-sized, appear to have fairly large extinc- 
tion cross-sections as well as single-scattering albedos at near- 
infrared wavelengths. At these wavelengths, they can thus con- 
tribute significantly to the scattering of radiation. The other types 
of particles appear to scatter less radiation. In particular the iron 
particle types have relatively small single-scattering albedos, 
while the silicate particle types have small albedos and small 
extinction cross-sections. For the silicate haze particles, the ex- 
tinction cross-section as calculated using the small-particle ap- 
proximation (see Section 2.1) also is plotted and can be seen to 
agree well with the results of the Mie-calculations, which con- 
vices us that our approximation of keeping the extinction cross- 
sections of these particles fixed and not only varying their optical 
thicknesses is valid. 



1500 K 




1 2 3 4 5 

Wavelength (micron) 



The product of the particle number density, its extinction 
cross-section at a given wavelength, and the geometric thick- 
ness of the atmospheric layer (assuming a homogeneous particle 
number density across the layer), yields the layer's extinction op- 
tical thickness at that wavelength. By multiplying this with the 
single-scattering albedo of the particles at the same wavelength, 
we obtain the layer's scattering optical thickness. The difference 
between the layer's extinction optical thickness and its scattering 
optical thickness is the layer's absorption optical thickness. 

2000 K 

10° ,:,:-. . -,:^ I 





10 








10 


c 







10 


t3 






10 


X 




LU 


10 




10 



1 2 3 4 5 

Wavelength (micron) 



o 1.0 I I 

"D 
0) 

S 0.8- 

D) /■ " " ~' - 

■i 0.6 -y' - 

g 0.4--.. - 

> ^ 
I, 0.2 - 

c _ ^ ^ 

1 2 ~ 3 4 5 

Wavelength (micron) 

Fig. 4. As Fig. [3] except for the 2000 K model atmosphere. 
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Fig. 3. Extinction cross-sections and single-scattering albedos 
of the four particle types for the 1500 K model atmosphere: 
the high silicate haze (solid lines); the low silicate haze (dot- 
ted lines); the iron cloud (dashed lines); and the AI2O3 cloud 
(dot-dashed lines). The thin solid line in the top panel shows the 
small-particle approximation for the high silicate haze, which 
agrees well with the results obtained using Mie calculations. 



Comparing Figs. [3] and |4] we see that the particles in the 
Tetf =2000 K atmosphere have lower albedos and are thus more 
absorbing than the particles in the reff=1500 K atmosphere. The 
reason is mostly that the former are smaller than the latter (see 
Fig.ID. 



4. Results 

Thermal emission spectra for the 1500 K case are shown in 
Fig. |5] The units plotted here are brightness temperature, since 
this does not depend on specific planet parameters like radius 
or distance to the observer Brightness temperature is the tem- 
perature corresponding to a black-body having the radiance that 
is observed or modeled. Brightness temperature has the addi- 
tional benefit that it roughly shows where in the atmosphere 
the radiance is coming from, since it can be compared with the 
real temperatures in the atmosphere. First a spectrum was cal- 
culated for an atmosphere with the drift-phoenix temperature 
profile and gas abundances, but without the clouds included. 
Although now the temperature profile is not self-consistent in 
the DRiFT-PHOENix framcwork, it will show how much including 
the clouds alone changes the emission spectrum, given a cer- 
tain temperature profile. The spectrum clearly shows the absorp- 
tion features of the water vapour, with wavelengths of low opac- 
ity giving higher brightness temperatures, because the thermal 
emission originates lower in the atmosphere. 

Adding in the scattering particles has a very strong effect on 
the spectrum. In the 1500 K case, the optical thickness of the 
particles is very large, with a total optical thickness of ~12 at 
a wavelength of 1 fim. As usual with small particles, the ex- 
tinction optical thickness decreases with increasing wavelength 
(see Fig. Is). Hence, the effect of the particles on the spectrum is 
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Fig. 5. Brightness temperature spectra for the 1500 K case. The 
solid line indicates the drift phoenix atmosphere without the par- 
ticles included. For the dotted line the scattering drift-phoenix 
particles are included. The dashed line shows the spectrum for 
purely absorbing particle, but same optical thickness (single- 
scattering albedo, or ssa in the legend, of zero). The dot-dashed 
line corresponds to the same particles, but now assumed per- 
fectly scattering (single-scattering albedo of unity). 



strongest at short wavelengths. The optical thickness of the parti- 
cles is larger than that of the water in the windows of low absorp- 
tion. If the particles are absorbing, this means that most of the 
radiation will come from the particles themselves at these wave- 
lengths, instead of the water lower in the atmosphere. Higher 
altitudes correspond here to lower temperatures, so the bright- 
ness temperature is reduced,as can be seen in Fig. |5] Because a 
more limited altitude range (and thus a more limited tempera- 
ture range in this case) is probed by the spectrum, the presence 
of the particles reduces the amplitude of the water features in the 
spectrum. 

We also calculated a spectrum with the single-scattering 
albedo of all particles set to zero and optical thicknesses are kept 
identical. This simulates the case when extinction is treated as 
pure absorption and scattering is neglected. The difference be- 
tween the dashed and dotted lines thus shows what the effect 
of the scattering is on the spectrum. In this case, the water fea- 
tures are slightly more subdued and the entire spectrum shifts 
by roughly 100 K in brightness temperature at low wavelengths 
compared to the scattering calculation. The effect of the scatter- 
ing is two-fold, corresponding to the two terms on the right-hand 
side of Eq. |3] Firstly, it allows more heat from the lower atmo- 
sphere to pass through than expected from extinction only. At 
the same time, the scattering particles locally reduce the emis- 
sivity of the atmosphere, leading to less radiation being emitted 
by that region of the atmosphere. Scattering can thus increase or 
reduce the emitted radiance compared to pure absorbing parti- 
cles, depending on whether the additional scattered radiation is 
more or less than the reduced emission. 

The spectra for the 2000 K case are shown in Fig.|6] Optical 
thicknesses of the cloud particles are lower for this case, which 
strongly reduces the effect of the particles on the spectrum. 
Particles also are generally more absorbing in this case, fur- 
ther reducing the effect of scattering. The difference between the 
spectrum that includes scattering and that without scattering is 
therefore small. 



Fig. 6. As Fig. |5] but now for the 2000 K case. 

In the above spectra the effect of scattering by the particles is 
significant, but not very large. This is because drift-phoenix pre- 
dicts small particles with a significant iron content throughout 
the upper atmosphere, making them relatively absorbing. For the 
1500 K case, particle extinction optical depth of unity is reached 
within the lower silicate haze layer (dotted lines in Fig. |3}, with 
significant extinction by the layers above this. Hence, the atmo- 
sphere that is emitting the thermal radiation is not very scat- 
tering in the drift-phoeni x model. However, p articles might be 
larger in other models (see lHelling et al. 2008ai , their Fig. 5), as 
well as in real atmospheres. Mie theory shows that larger parti- 
cles of the same composition will generally have higher single- 
scattering albedos and more forward-scattering phase functions 
than small particles. This effect of the particle size also is visi- 
ble in the single-scattering albedos in Figs. |3]and|4] The more 
forward-scattering particles with higher single-scattering albe- 
dos would result in more radiation from lower altitudes being 
scattered upward into space. Furthermore, drift-phoenix has a 
significant iron content in its particles in the upper atmosphere, 
which is strongly absorbing. If the iron content is lowered, the 
single scattering albedo also quickly rises. To illustrate the effect 
of highly scattering particles, we also calculated a case where 
the particle single-scattering albedos are set to unity (i.e. non- 
absorbing particles). This is the extreme case, but it is realistic 
for particles larger than the wavelength with a low iron content. 
For ease of comparison between the cases with different single- 
scattering albedos, the phase functions and optical thicknesses 
are kept identical here. Figs.|5]and|6]show this case as well (dot- 
dashed lines). It is clear that the scattering makes a very large dif- 
ference in how the spectrum looks here: gas absorption features 
are much more apparent and overal brightness temperature levels 
are changed substantially. Also note that the relative strengths of 
the water bands can change significantly by changing the single- 
scattering albedo of the particles. Thus, if clouds would not be 
considered in the interpretation of this thermal emission spec- 
trum, one could derive completely different temperature profiles 
and gas adundances than what is actually present in the atmo- 
sphere. If phase functions also would be more forward scatter- 
ing, as is generally the case with larger particles, this difference 
with non-scattering particles will be even more pronounced. 

To illustrate the potential effect of a more forward-scattering 
phase function, we calculate two spectra with only a difference 
in the phase function. The temperature and optical thicknesses 
at 1 micron are taken from the 1500 K case, but instead of 
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the particle wavelength-dependence of Fig. |3] we assume optical 
thickness to be constant with wavelength and single-scattering 
albedos of unity, as is common for large non-absorbing parti- 
cles. Also, all particles throughout the atmosphere were assumed 
identical. The two different phase functions were also assumed 
constant with wavelength and were taken from the Teff = 1500 
K AI2O3 cloud at a wavelength of 0.3 micron (asymmetry pa- 
rameter of g = 0.83) and the Ten =2000 K high silicate haze 
at 5 micron = 3 • 10 The first is thus strongly forward- 
scattering, whereas the latter is almost isotropically scattering. 
The two spectra with the two different particles are shown in 
Fig. I2] Even though the particle optical thicknesses for both 
cases are identical, the spectra differ very strongly. The forward- 
scattering particles allow more heat from the warmer region s be- 
low the clouds to escape to space, giving higher radiances and 
stronger absorption features. Hence, assuming isotropic scatter- 
ing can give errors of many tens of percent on the spectrum if 
the particles are in fact strongly forward-scattering. 
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Fig. 7. Spectra for two atmospheres with temperatures and op- 
tical thicknesses as the 1500 K case at 1 fim (now taken to be 
constant with wavelength) and single-scattering albedos of unity, 
with a difference in phase function. The solid line is the spectrum 
for the more forward-scattering particles {g = 0.83), the dashed 
line is for the isotropically scattering particles (g - 3 ■ 10""*). 



5. Conclusions and discussion 

Cloud and/or haze particles in the atmospheres of hot exoplan- 
ets or brown dwarfs can have a strong effect on their thermal 
emission spectra, changing their brightness and colour at differ- 
ent wavelengths. These thermal emission spectra can be mea- 
sured either directly, or from their secondary eclipse. The effect 
of the clouds is strongest at shorter wavelengths and colder tem- 
peratures for the DRIFT-PHOENIX model, and generally for small 
particles. This makes the brown dwarf more 'red' in the in- 
frared. Although the latter effect is clear from previous studies, 
not much attention has been given in the past to the contribu- 
tion from scattered light to the spectrum. Here, we show that 
not only the extinction of the cloud, but also the scattering na- 
ture of the particles can be important in determining the emis- 
sion spectra for these objects. Scattering can affect the emission 
spectra especially when the particles are large and have little 
iron in them. Hence, knowledge of the scattering properties such 



as single-scattering albedo and phase function can be crucial in 
calculating accurate spectra of sub-stellar atmospheres. Errors 
in these parameters will result in errors in atmospheric prop- 
erties derived from infrared direct measurements or secondary 
eclipses. Also the self-consistent calculation of the temperature 
structure in such an atmosphere can thus depend significantly 
on the assumed or calculated scattering properties in the atmo- 
sphere. For instance. Fig. |7] shows that more forward-scattering 
phase functions allow more heat from lower altitudes to escape 
to space, leading to an increase in cooling of these lower re- 
gions compared to more isotropically scattering particles. This 
also will affect calculations made using the isotropically scat- 
tering PHOENIX model, if the particles are more scattering than 
the silicate haze in the two drift-phoenix cases presented here. 
The errors due to the parameterisation of the phase function in 
terms of g are bound to be less severe, but probably can reach 
several percent in the calculated radiance. Thus, care must be 
taken in the choice of radiative transfer method for heat balance 
or spectral calculations, depending on the type of particles that 
are predicted. 

Current analyses of transiting exoplanet emission spectra, as 
derived from secondary eclipses, suffer from a large degeneracy, 
with a wide range of atm ospheres fitting the limited me asure- 
ments equally well (e.g. iMadhusudhan & Seager 20091) . This 
means that atmospheric properties like temperature and compo- 
sition are poorly constrained by the observations. Most current 
secondary eclipse measurements are taken in the near-infrared 
by the warm Spitzer mission and ground-based telescopes and 
none of the current fits to the measurements include clouds. We 
show here that adding scattering clouds only will add to the 
degeneracy, as e.g. the depths of absorption features or bright- 
ness temperatures at a certain wavelength in the near-infrared 
strongly depend on the scattering properties and optical thick- 
nesses of the particles, as well as the temperature structure and 
gas abundances. As there is no evidence that the currently anal- 
ysed transiting exoplanets are cloud-free, clouds should be taken 
into account when fitting the observations, hence increasing the 
uncertainty in the other free parameters such as temperature 
and gas abundances. A cloud physics model like drift-phoenix 
will be valuable in exploring the parameter space and physically 
explaining the appearance of the measured spectra. Also, the 
best-fit temperature profiles and gas adundances, derived with- 
out clouds, can be substantially in error when clouds are present 
in the atmosphere. 

Our spectra also show that the composition of clouds them- 
selves cannot be constrained well from near-infrared measure- 
ments of the thermal emission alone, since multiple combina- 
tions of composition mixtures and particle sizes can give very 
similar extinction and scattering properties. Measurements at 
longer wavelengths, where the expected silicate particles have 
absorption features, can possibly help to reduce the degener- 
acy, for both directly dete cted exoplanets and transiting ex- 
oplanets (He lUng et al. 2008c). Also observations of reflected 
starligh t at short wavelengths and transit transmission mea sure- 
ments dLecavelier Pes Etangs et al. 20081: iPont et al. 20081) can 
help to determine the haze and cloud structure and so more ac- 
curately determine the temperature structure and gas adundances 
from secondary eclipse observations. Note that scattering of the 
starlight may play a role in transit transmission measurements as 
well. 

The fact that the thermal radiation can be scattered 
also means that it can be polarised. For a symmetri- 
cal, spherical planet the disc-integrated polarisation will be 
zero, but when the planet or brown dwarf is flattened 



7 



de Kok et al.: Scattering of thermal radiation on brown dwarfs and hot exoplanets 



dSengupta & Marlev 200911 20101) or horizontally inhomoge- 
neous it can result in the infrared thermal radiation being po- 
larised. Our radiative transfer code also includes an accurate de- 
scription of the polarisation of the light and so it can calculate 
polarisation of the light given a certain viewing geometry. In a 
separate work we will investigate the possible polarisation sig- 
nals from the thermal radiation of horizontally inhomogeneous 
planets. 

Finally, it also is interesting to note that the small dark 
particles predicted by drift-phoenix also will result in a lower 
geometric albedo of the planet at visible wavel engths com- 
pared to Rayleigh scattering (Burrows et al. 2008b) or relativel y 
large enstatite particles (Sudarskv et al. 2000; Hood et al. 20081) . 
The latter two do not model the cloud physics, but assume 
an ad hoc par ticle size distribution. The radiative transfer in 
iHood et al. (2 008) differs from other work cited here, as it 
takes into account all three spatial dimensions of the atmo- 
sphere. Because the small iron-rich particles are dark, this haze 
may contribute t o the low albedos of some hot Jupiters (e.g. 
iRowe et al. 20081) . The darkening effect of small iron particles 
has already been noted by Seager et al. (2000), but the drift- 
phoenix model might provide a physical basis for this argument. 
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